SU-4240-601 
|hep-ph/95014T7 
January 31, 1995 



Exploring tttt scattering in the — picture 



c 



Francesco Sannino aib,c f\ and Joseph Schechter a Q 



a Department of Physics, Syracuse University, Syracuse, New York, 132^-1130. 



b 



Dipartimento di Scienze Fisiche, Mostra d'Oltremare Pad. 19, 80125 Napoli, Ita- 
lia. 

Istituto Nazionale di Fisica Nucleare, Sezione di Napoli, Mostra d'Oltremare, 
Pad. 19 1-80125 Napoli, Italia. 

Abstract 

In the large N c approximation to QCD, the leading 7T7t scattering amplitude is ex- 
pressed as the sum of an infinite number of tree diagrams. We investigate the possibility 
that an adequate approximation at energies up to somewhat more than one GeV can be 
made by keeping diagrams which involve the exchange of resonances in this energy range 
in addition to the simplest chiral contact terms. In this approach crossing symmetry is 
automatic but individual terms tend to drastically violate partial wave unitarity. We first 
note that the introduction of the p meson in a chirally invariant manner substantially de- 
lays the onset of drastic unitarity violation which would be present for the current algebra 
term alone. This suggests a possibility of local (in energy) cancellation which we then 
explore in a phenomenological way. We include exchanges of leading resonances up to the 
1.3 GeV region. However, unitarity requires more structure which we model by a four 
derivative contact term or by a low lying scalar resonance which is presumably subleading 

in the — expansion, but may nevertheless be important. The latter two flavor model 

^ c 

gives a reasonable description of the phase shift 5® up until around 860 MeV, before the 
effects associated which the KK threshold come into play. 
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1 Introduction 



The study of tttt scattering has been one of the classical methods for investigating the nature 
of the strong interaction. Many elegant ideas have been proposed At the present time the 
standard approach at low energies is based on chiral perturbation theory ||. This enables 
one to nicely understand the scattering amplitudes near threshold(< 400 MeV). However 
it is very difficult to extend this treatment to higher energies since the pole structure of the 
resonances in this region cannot be easily reproduced by a truncated power series expansion 
in energy. Our interest in this paper will be to investigate in a schematic sense how the 
description of tttt scattering might be extended up to slightly past the 1 GeV region in a new 
chiral picture. 

It has been clear for many years that in the region just beyond the threshold region the 
effects of p exchange dominate. However in the chiral perturbation program the effects of 
the p arise in the second order of the energy expansion @ . This is of course due to the fact 
that the usual chiral program is mainly devoted to improving the description of the dynamics 
in the threshold region in which the p does not explicitly appear. For going beyond the 
threshold region we would like an approach which can treat the p and other resonances at 
the first stage of an iterative procedure. 

Such an approach is suggested by the large iV c approximation to QCD. As reviewed in [|j 
for example the leading order — — approximation to tttt scattering is obtained by summing all 
possible tree diagrams corresponding to some effective Lagrangian which includes an infinite 
number j|] of bosonic resonances of each possible spin. In addition it is allowed to include 
all possible contact terms. This clearly has the right structure but initially seems to be so 
general as to be practically useless. Here we will argue that this may not be the case. 
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An amplitude constructed according to the above prescription will automatically satisfy 
crossing symmetry. On the other hand just calculating the tree approximation to an effective 
Lagrangian will not guarantee that unitarity is satisfied. This is the handle we will use 
to try to investigate additional structure. Unitarity has of course the consequence that the 
amplitude must have some suitable imaginary term which in the usual field theory is provided 
by loop diagrams. However the leading — approximation will give a purely real amplitude 
away from the singularities at the direct s — channel poles. We may consider the imaginary 
part of the leading — amplitude to consist just of the sum of delta functions at each such 
singularity. Clearly, the real part has a much more interesting structure and we will mainly 
confine out attention to it. Furthermore we will assume that the singularities in the real part 
are regularized in a conventional way. 

Unitarity has the further consequence that the real parts of the partial wave amplitudes 
must satisfy certain well known bounds. The crucial question is how these bounds get satisfied 
since, as we will see, individual contributions tend to violate them badly. At first one might 
expect that all of the infinite number of resonances are really needed to obtain cancellations. 
However the success of chiral dynamics at very low energies where none of the resonances have 
been taken into account suggests that this might not be the case. At the very lowest energy the 
theory is described by a chiral invariant contact interaction which, however, quite soon badly 
violates the unitarity bound. It will be observed that the p exchange tames this bad behavior 
dramatically so that the unitarity bound is not badly broken until around an energy beyond 
2 GeV. This suggests that there is a local cancellation between resonances which enforces the 
bound. The local cancellation is not easily predicted but if true in general it would greatly 
simplify the task of extending the phenomenological description of scattering processes to 
higher energies. Including just the effects of the p and the ir particles corresponds to including 
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just the s — wave quark antiquark states in our model. A natural next step , which we shall 
also explore here, would be to include the exchange of the allowed p — wave quark antiquark 
states. At the same time the quark model suggests that we include the first radial excitations 
of the s-wave states which in fact lie near the p — wave states. Theoretically glueball and 
exotic states are suppressed in irir scattering according to the large N c approximation 

In our analysis the pions will be treated as approximate Goldstone bosons corresponding 
to the assumption that the theory has a spontaneously broken chiral symmetry. Actually 
the need for spontaneous breakdown of this symmetry can be argued from the — approach 
itself ||. The effective Lagrangian we use shall be constructed to respect this symmetry. 
Furthermore for the purpose of the initial exploration being performed here we shall consider 
resonance interaction terms with the minimum number of derivatives and shall also neglect 
chiral symmetry breaking terms involving the resonances. When going beyond the initial 
stage we will be forced to proceed in a more phenomenological way. 

In section 2, after the presentation of the partial wave amplitudes of interest, we show how 
the introduction of the p— meson in a chirally invariant manner substantially delays the onset 
of the severe unitarity violation which would be present in the simplest chiral lagrangian of 
pions. The program suggested by this local cancellation is sketched. 

Section 3 is concerned with the contribution to the pion scattering amplitude from the 
next group of resonances - those in the range of the p — wave qq bound states in the quark 
model. It is observed that a four derivative contact term can be used to restore unitarity up 
to about 1 GeV. 

In section 4, a possibly more physical way to restore unitarity is presented which makes 
use of a — subleading contribution due to a very low mass scalar (presumably qqqq) state. 
Finally, section 5 contains a brief summary and discussion of some directions for future 
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work. 



2 Current algebra plus p exchange 

In this section we will study the partial waves for tttt scattering computed in a chiral La- 
grangian model which contains both the pseudoscalar and vector mesons, (i.e., the lowest 
lying s-wave quark antiquark bound states). 

The kinematics are discussed in Appendix A, where the partial wave amplitudes Tf are 
defined. They have the convenient decomposition: 



where 8j (s) are the phase shifts and 77/ (s) (satisfying < 77/ (s) < 1) are the elasticity 
parameters. Extracting the real and imaginary parts via 






rjf sin(25j) 
2 ' 
1-rif cos(25f) 
2 




(2.3) 



leads to the very important bounds 




< If < 1. 



(2.4) 



For fixed r\\ the real and imaginary parts lie on the well known circle in the Argand-plane 

1 T 2 

R,{ 2 + (// — -) 2 = ■ This formula also enables us to solve for l{ as : 




(2.5) 



Let us use (^J5|) for an initial orientation. Near threshold r\\ = 1, R{ is small and we should 





(2.6) 
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In the large N c limit the amplitude near threshold is purely real and of the order — . This 

1 

is consistent with ( |2.6| ) which shows that l( (s) is of order an d hence comes in at the 
second order. This agrees with the chiral perturbation theory approach 0] in which Rj(s) 
comes from the lowest order tree diagram while if arises from the next order loop diagram. 
On the other hand, when we depart from the threshold region the — approach treats the 
contribution of the p-meson at first order while the chiral perturbation theory approach treats 
it at second and higher orders. 

Now, pion physics at very low energies is described by the effective chiral Lagrangian, 

L x = -^Tr(d^Ud^) + Tr(B(U + t/t)) (2.7) 
o 

wherein U = e F ^ and <f> is the 3x3 matrix of pseudoscalar fields. F n = 132 MeV is 

F 2 m 2 

the pion decay constant. Furthermore B = diag(Bi, B\, B%) , where B\ = 71 ^ 71 and 

B3 = — , describes the minimal symmetry breaking. We shall choose = 

137 MeV. 

A straightforward computation using (2/7) yields the it-it scattering amplitude @ defined 



in (A.l): 



(s — mV 



A(s,t,u)=2 y p2 (2i 



This equation will be called the current algebra result. With (A.2)-(A.4) we obtain -R{j(s) = 
Tq(s) as illustrated in Fig. 1. The experimental Roy curves || are also shown. Up till 
about 0.5 GeV the agreement is quite reasonable (and can be fine tuned with second order 
chiral perturbation terms) but beyond this point Rq keeps increasing monotonically and badly 
violates the unitarity bound (2.4). We will see that the introduction of the p-meson greatly 



improves the situation. 
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There are several different but essentially equivalent ways to introduce vector mesons into 
the chiral invariant Lagrangian. A simple way Q is to treat the vectors as gauge particles 
of the chiral group and then break the local symmetry by introducing mass-type terms. The 
3x3 matrix of the vector fields, p^ is related to auxiliary linearly transforming gauge fields 

and by 

Ajl = tPvtf + Ud^ (2.9) 
>> g 

A* = &PnZ + U%£, (2-10) 

i 

where £ = U 2 and g is a dimensionless coupling constant. Under a chiral transformation 
4l0> 



u u L uul 



£ - U L £tf = K£U R (2.11 
(which also defines the matrix K(cj>, Ul, Ur) and p^ behaves as 

Pf 1 



It is convenient to define 



Kp^rf + -Kd^rf. (2.12) 
9 



v» = l(tdv? + t%t) (2-13) 
Pa. = (2-14) 



Vll -> Ap^jrt (2.i5) 

^ -» Kv^ + iK3^. (2.16) 

These quantities enable us to easily construct chiral invariants and will also be useful later 
[11|. The chiral Lagrangian including both pseudoscalars and vectors that one gets can be 



which transform as 



rewritten as the sum of L\, in Q2.7| ) and the following: 
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1 m n r 

L 2 = —Tr(Fp U {p)Fp V {p)) - ^Tr [(g PfJ 



(2.17) 



where Fp U (p) = dpp u — d u p p — ig[pp, p v \- The coupling constant g is related to the p- meson 
width by 



T(p -> 2tt) 



12iriJip 



m. 



9p-n 



(2.18) 



We choose m p = 0.769 GeV and g p7T7T = 8.56. Symmetry breaking contributions involving 
the p are given elsewhere [|l2| but are small and will be neglected here. The Lagrangian piece 



in (2.17) yields both a pole-type contribution (from the PpVp cross term) and a contact term 
contribution (from the VpVp term) to the amplitude at tree level ||]: 



A(s,t,u) = (2. 



2 \ml-t 



s t — s 
+ ^ 



— u 



3 pirn 
2777,2 



[(t -s) + (u- s)\ 



(2.19) 



We notice that the entire second-term in ( |2.17| ) is chiral invariant since v p and gpp transform 
identically. However the Tr(ppVp) and Tr(vpVp) pieces are not separately chiral invariant. 
This shows that the addition of the p meson in a chiral invariant manner necessarily introduces 
a contact term in addition to the minimal pole term. Adding up the terms in (2.19) yields 
finally 



A(s, t, u) 



m: 



F 2 



Q 2 

2m 2 



t(u — s) u(t — s) 

o 7~ o 



TTlf 



777^ — U 



(2.20) 



In this form we see that the threshold (current algebra) results are unaffected since the 
second term drops out at t = u = 0. An alternative approach [|13| to obtaining ( 2.20| ) 
involves introducing a chiral invariant pinr interaction with two more derivatives. 

A(s,t,u) has no singularities in the physical region. Reference to (A. 2) shows that the 
isospin amplitudes T° and T 2 also have no singularities. However the T 1 amplitude has 
the expected singularity at s = mi. This may be cured in a conventional way, while still 
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maintaining crossing symmetry, by the replacements 

1 1 



t, u m 2 — t,u — irripT 



(2.21) 



in ( |2.20| ) 0. A modification of this sort would enter automatically if we were to carry the com- 
putation to order — However we shall regard ( [2.21D as a phenomenological regularization 
of the leading amplitude. 

Now let us look at the actual behavior of the real parts of the partial wave amplitudes. 
Rq, as obtained from ( p.20|) with ( |2 . 2 1| ) , is graphed in Fig. 2 for an extensive range of \/s, 
together with the pions only result from fl2,8| ). We immediately see that there is a remarkable 
improvement; the effect of adding p is to bend back the rising i?§( s ) so there is no longer a 
drastic violation of the unitarity bound until after y/s = 2 GeV. There is still a relatively 
small violation which we will discuss later. Note that the modification ( [2.21 ) plays no role in 
the improvement since it is only the non singular t and u channel exchange diagrams which 
contribute. 

It is easy to see that the delayed drastic violation of the unitarity bound \Rj\ < — is 
a property of all partial waves. We have already learned from ( |2.20| ) that the amplitude 
A(s, t, u) starts out rising linearly with s. Now fl2.19|) and ( |A.3| ) show that for large s the p 
exchange terms behave as s°. The leading large s behavior will therefore come from the sum 
of the original current- algebra term and the new contact-term: 

2s 3k\ , m 2 p 



But k is numerically around 2 H], so A(s,t,u) eventually decreases linearly with s. This 
turn-around, which is due to the contact term that enforces chiral symmetry, delays the onset 
of drastic unitarity violation until well after the p mass. It thus seems natural to speculate 



One gets a slightly different results if the the regularization is applied to (2.19). 



that, as we go up in energy, the leading tree contributions from the resonances we encounter 
( including both crossed channel as well as s-channel exchange) conspire to keep the Rj(s) 
within the unitarity bound. We will call this possibility, which would require that additional 
resonances beyond the p come into play when R{(s) from (2.20) start getting out of bound, 
local cancellation. 



In Fig. 3 we show the partial waves R\ and l\ computed using and (fl2~ll ). Not 



surprisingly, these display the standard resonant forms. For completeness we present the Rq 
and i?2 amplitudes in Fig. 4. We may summarize by saying that the results of this section 
suggest investigating the following recipe for a reasonable approximation to the tttt scattering 
amplitude up to a certain scale E max . 

1. Include all resonances whose masses are less than E max + A, where A ~ several-hundred 
MeV. This express the hoped for local cancellation property. 

2. Construct all possible chiral invariants which can contribute, presumably using the 
minimal number of derivatives. Compute all tttt —* tttt tree diagrams, including contact 
terms. Regulate the resonance denominators in a manner similar to ( |2.21| ), but restrict 
attention to the real part. Interpret the manifestly crossing symmetric result as the 
leading order in — real tttt scattering amplitude. 

3. Obtain the imaginary parts of the partial wave amplitudes, using fl2.5|) . The rjf (s) 
might be computed by including channels other than tttt. 

We will start to explore this program by checking whether the inclusion of the next group 
of resonances does enable us to satisfy the unitary bound for Rq. For simplicity we restrict 
ourselves to a two-flavor framework. It is straightforward to generalize the scheme to three 
flavors. 



3 The next group of resonances 

To keep our investigation manageable we shall mainly restrict attention to the partial wave 
amplitude Rq(s). As we saw in the last section, this is the one most likely to violate the 
unitarity bound. The first task is to find the effective lagrangian which should be added to 
(2.7) plus ( p. 17 ). There is no a priori reason not to add chiral invariant tttt contact interactions 



with more than two derivatives. But the most characteristic feature, of course, is the pionic 

interactions of the new resonances in the energy range of interest, here up to somewhat more 

1 

than 1.0 GeV. Which resonances should be included ? In the leading — approximation, qq 
mesons are H ideally mixed nonets, assuming three light quark flavors. Furthermore, the 



exchange of glueballs and exotic mesons are suppressed in interactions with the qq mesons. 
The — approximation thus directs our attention to the p — wave qq resonances as well as 
the radial excitations of the s — wave qq resonances. 



The neutral members of the p — wave qq nonets have the quantum numbers J = ++ , 
1 ++ , 1 H and 2 ++ . Of course, the neutral members of the radially excited s — wave qq nonets 
have J PC = ^ and 1 . Only members of the ++ , 1 and 2 ++ nonets can couple to 
two pseudoscalars Q. By G-parity conservation we finally note that it is the 1 = member of 
the ++ and 2 ++ nonets and the 1=1 member of the 1 nonet which can couple to two 
pions. Are there good experimental candidates for these three particles ? 

The cleanest case is the lighter 1 = member of the 2 ++ nonet; the /2(1270) has, 
according to the August 1994 Review of Particle Properties (RPP) [|l^], the right quantum 
numbers, a mass of 1275 ± 5 MeV, a width of 185 db 20 MeV, a branching ratio of 85% 

into two pions, and a branching ratio of only 5% into KK. On the other hand the ^2(1525) 

4 It is possible to write down a two point mixing interaction between v and radially excited 0~ + particles 
etc., but we shall neglect such effects here. 
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has a 1% branching ratio into tttt and a 71% branching ratio into KK. It seems reasonable 
to approximate the 2 ++ nonet as an ideally mixed one and to regard the /2(1270) as its 
non-strange member. 

The /o(1450) is the lightest listed |l^] particle which is a candidate for a radial excitation 
of the usual p(770). It has a less than 1% branching ratio into KK but the tttt branching 
ratio, while presumably dominant, is not yet known. With this caution, we shall use the 
p(1450). The p(1700) is a little too high for our region of interest. 

An understanding of the 1 = 0, ++ channel has been elusive despite much work. The 
RPP flTjl gives two low lying candidates: the /o(980) which has a 22% branching ratio into 
KK even though its central mass is below the KK threshold and the /o(1300) which has about 
a 93% branching ratio into tttt and a 7% branching ratio into KK. We shall use the /o(1300) 
here. It is hard to understand why, if the /o(980) is the ss member of a conventional ++ 
nonet, it is lighter than the /o(1300). Most likely, the /o(980) is an exotic or a KK molecule 
]l6[ . If that is the case, its coupling to two pions ought to be suppressed in the — picture. 
This is experimentally not necessarily true but we will postpone a discussion of the /o(980) 
as well as other possible light ++ resonances to the next section. 

Now we will give, in turn, the tttt scattering amplitudes due to the exchange of the 
/ (1300), the /2(1270) and the p(1450). 

3.1 The / (1300) 

Denoting a 3 x 3 matrix of scalar fields by S we require that it transform as S — > KSK^, 
(see ( |2.1lD ) under, the chiral group. A suitable chiral invariant interaction, using ( 2.14 ), is 



L f0 = - lo F%Tr(S P ^) = -joTriSd^d^) + ■■■ (3.1) 
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where the expansion of £ was used in the second step. It is interesting to note that, in the 
present formalism, chiral symmetry demands that the minimal S(f)(f) interaction must have 
two derivatives. Specializing to the particles of interest and taking /q to be ideally mixed 
leads to 

L f0 = -^(d^-d^) + .... (3.2) 
The partial width for /q(1300) — * tttt is then 
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roMiMO) _ m) _ (, _ M) 5 (M% - 2,^ ,3.3) 

The RPPg|| lists r toi (/ (1300)) = 0.15 - 0.40 GeV and M fo = 1.0 - 1.5 GeV. For def- 
initeness we shall choose r tot (/o(1300)) = 0.275 GeV and Mf = 1.3 GeV. These yield 
|7o| = 2.88 GeV -1 . Using ( |3.2| ) we find the contribution of /o exchange to the tttt scattering 
amplitude, defined in ( |A.1| ), to be: 

A„Xs, t ,u) = f^0. (3.4) 

JO 



Actually, as discussed around ( |2.2l|) , the singularity in the real part of (p.4[) will be regulated 
by the replacement 

1 M l - a 
- & (3.5) 



M l- S (Ml-sY + Mlp- 
3.2 The / 2 (1270) 

We represent the 3x3 matrix of tensor fields by T^ v (satisfying T^ v = T Vfi , and = 0) 
which is taken to behave as T^ v — ► KT^ V K^ under chiral transformation. A suitable chiral 
invariant interaction is 

L T = -72^ 2 Tr(l> W ). (3.6) 
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Specializing to the particles of interest, this becomes 

L 



-|(/2V(^7f-^7f) + ---. 



(3.7) 



In this case we note that the chiral invariant interaction is just the same as the minimal one 
we would have written down without using chiral symmetry. The partial width is then 



r(/ 2 (1270) -» vrvr) 



72 Pi 
20vr Mf 



(3.8) 



where is the pion momentum in the fi rest frame. This leads to I72I = 13.1 GeV 



-1 



To calculate the /2 exchange diagram we need the spin 2 propagator |17|| 

"1 



M%+q 2 



1 

3 



where 



Kv + 



M 2 
h 



(3.9) 



(3.10) 



A straightforward computation then yields the j% contribution to the tttt scattering amplitude: 



Af 2 (s,t,u) 



7l 



2(M 2 -s) \ Y m * + Y m - S -3 S + 2 {t +U ' 



2 m\s 2 

3 M 2 

h 



s 3 s 4 
+ 



6M 2 6Mj ' 

h J2 / 



(3.11) 



Again the singularity will be regulated as in (3.5). 
3.3 The p(1450) 

We may read off the proper chiral invariant contribution^ of the p(1450) to the tttt scattering 
amplitude from the second term of ( |2.20| ) 



A p '(s, t, u) 



2m 2 , 



t(u — s) u(t — s) 

2 T ' 2 

m p , — t m p , — u 



(3.12) 



5 It is not necessary to introduce the p(1450) as a massive gauge field as we did for the p(770), but the 
answer is the same. See [Ol for further discussions. 



13 



where g p i wn is related to the p' — > tttt partial width by 

We shall use (^'^1 = 7.9 corresponding to T(p — > 27r) = 288 MeV. 
3.4 / (1300)+/ a (1270)+p(1450) 

Now we are in a position to appraise the contribution to R® of the next group of resonances. 
This is obtained by adding up (jOD, fmD and fTT^ ) and using flA^) - (gjg) . The individual 



pieces are shown in Fig. 5. 

Note that the /o(1300) piece is not the largest, as one might at first expect. That honor 
goes to the contribution which is shown divided into the s-channel pole piece and the 
(t + u) pole piece. We observe that the s-channel pole piece, associated with the /2, vanishes 
at y/s = Mf 2 . This happens because the numerator of the propagator in ( |Q| ) is precisely a 
spin 2 projection operator at that point. The p(1450) contribution is solely due to the t and 
u channel poles. It tends to cancel the t and u channel pole contributions of the /2(1270) 
but does not quite succeed. The t and u channel pole contributions of the /o(1300) turn out 
to be negligible. Notice the difference in characteristic shapes of the s and (t + u) exchange 
curves. Fig. 6 shows the sum of all these individual contributions. There does seem to be 
cancellation. At the high end, Rq starts to run negative well past the unitarity bound fl2.4| ) 
around 1.5 GeV. But it is reasonable to expect resonances in the 1.5 — 2.0 GeV region to 
modify this. The maximum positive value of Rq is about 1 at yfs = 1.2 GeV . This would 
be acceptable if the tt + p contribution displayed in Fig. 2, and which must be added to the 
curve of Fig. 6, were somewhat negative at this point. However this is seen not to be the 
case, so some extra ingredient is required. The — approach still allows us the freedom of 
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adding four, and higher derivative contact terms. More physically, there is known to be a 
rather non trivial structure below 1 GeV in the I = J = channel. 

3.5 4— Derivative Contact Terms. 

First let us experiment with four-derivative contact terms. So far we have not introduced 
any arbitrary parameters but now we will be forced to do so. There are two four-derivative 
chiral invariant contact interactions which are single traces in flavor space: 

L 4 = a Tr{d^Ud v U ] d^Ud u U ] ) + b Tr^Ud^ d v Ud u U ] ) (3.14) 

where a and b are real constants. The single traces should be leading in the — expansion. 
Notice that the magnitudes of a and b will differ from those in the chiral perturbation theory 
approach j|] since the latter essentially also include the effects of expanding the p exchange 
amplitude up to order s 2 . The four pion terms which result from (3.14) are: 



2a (dfj/K ■ d^) 2 + (6 - o) {d„7f ■ d^f] +■■■. (3.15) 



This leads to the contribution to the -kit amplitude: 

A 4 (s, t,u) = ^ \a ((t - 2ml) 2 + (« " ^l) 2 ) + (b - a)(s - 2ml) 2 } ■ (3-16) 



i?4 

± 7T 

Plausibly, but somewhat arbitrarily, we will require that ( 3.16| ) yields no correction at thresh- 



old, i.e. at s = 4m^, t = u = 0. This gives the condition b = —a and leaves the single param- 



eter a to play with. In Fig. 7 we show Rq, as gotten by adding the piece obtained from ( 3.16 ) 



for several values of a to the contribution of ir + p, plus that of the next group of resonances. 
For a = +1.0 x 10 -3 the four-derivative contact term can pull the curve for down to avoid 
violation of the unitarity bound until around = 1.0 GeV. The price to be paid is that 
Rq decreases very rapidly beyond this point. We consider this to be an undesirable feature 
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since it would make a possible local cancellation scheme very unstable. Another drawback of 
the four-derivative contact term scheme is that it lowers Rq(s) just above threshold, taking it 
further away from the Roy curves. Let us therefore set aside the possibility of four and higher 
derivative contact terms and try to find a solution to the problem of keeping Rq within the 
unitarity bounds in a different and more phenomenological way. 

4 Low energy structure 



Let us investigate the addition of low energy [18| exotic states whose contributions to tttt 
scattering should be formally suppressed in the large N c limit. Experimentally we know that 
there is at least one candidate - the /o(980) mentioned in the previous section. According to 



the RPP [15] its width is 40 — AOOMeV and its branching ratio into two pions is 78.1 ± 2.4%. 
To get an idea of its effect we will choose T to t = ^OMeV and use the formulas fl3.2|) -([3li|) 
with the appropriate parameters. In Fig. 8 we show the sum of the contributions of the next 
group added to tt + p with the effect of including the /o(980). It is seen that the /o(980) does 
not help unitarity - below y/s = 980MeV it makes the situation a little worse while above it 
improves the picture slightly. 

What is needed to restore unitarity over the full range of interest and to give better 
agreement with the experimental data for yfs < 900 MeV? 

i. Below 450 MeV, Rq(s) actually lies a little below the Roy curves. Hence it would be 
nice to find a tree level mechanism which yields a small positive addition in this region. 

tt. In the 600 — 1300 MeV range, an increasingly negative contribution is clearly required 
to keep Rq within the unitarity bound. 

It is possible to satisfy both of these criteria by introducing a broad scalar resonance (like 
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the old a) with a mass around 530 MeV. Its contribution to A(s, t, u) would be of the form 
shown in ([O]) and ( |3.5| ) which we may rewrite as: 

32vr G {s-2ml) 2 (Ml- s) 



3H Ml (s - Ml) 2 + MIG' 



2 ' 



(4.1) 



7 2 327r G ( m 2 \^ ( m 2 \ 2 

where we have set y —-^3, r -> G' and H = I 1 - 4-^| j I 1 - 2-^| I is approx- 
imately one. If this were a typical resonance which was narrow compared to its mass and 
which completely dominated the amplitude, we would set G ~ G' « V. However for a very 
broad resonance it may be reasonable to regard G' as a phenomenological parameter which 

could be considered as a regulator in the sense we have been using. Choosing M a = 0.53 GeV, 
G 

— = 0.31 and G 1 = 380 MeV, the contribution to Rq of ( |4.1| ) is shown in Fig. 9. The curve 
G 

goes through zero near 0.53 GeV (there is a small shift due to the crossed terms). Below 
this value of y/s it adds slightly in accordance with point i while above 0.53 GeV it sub- 
tracts substantially in the manner required by point ii. This is the motivation behind our 
choice of M a = 0.53 GeV . Adding everything - namely the it + p piece, the next group piece 



together with the contributions from (4.1) and the /o(980) - results in the curves shown in 
Fig. 10 for three values of G' . These curves for Rq(s) satisfy the unitarity bound \Rq\ < ^ 
until y/s w 1.3 GeV. After 1.3 GeV, the curves drop less precipitously than those for the 
four-derivative contact term in Fig. 7. 

Fig. 10 demonstrates that the proposed local cancellation of the various resonance ex- 
change terms is in fact possible as a means of maintaining the unitarity bound for the (by 
construction) crossing symmetric real part of the tree amplitude. Essentially, just the three 
parameters M a , G and G' have been varied to obtain this. The other parameters were all 
taken from experiment; when there were large experimental uncertainties, we just selected 
typical values and made no attempt to fine-tune. Procedurally, G', G and M a were adjusted 
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to obtain a best fit to the Roy Curves below 700 MeV; this turned out to be what was needed 

for unitarity beyond 700 MeV. It was found that M a had to lie in the 530 ± 30 MeV range 
G' 

and that — had to be in the 0.31 ± 0.06 range in order to achieve a fit. On the other hand 
G' could be varied in the larger range 500 ± 315 MeV. It is also interesting to notice that 
the main effect of the sigma particle comes from its tail in Fig. 9. Near the pole region, its 
effect is hidden by the dominant ir + p contribution. This provides a possible explanation of 
why such a state may have escaped definitive identification. For the purpose of comparison 
we show in Fig. 11, the total Rq together with the 7r + p and current- algebra curves in the 
low energy region. 

It is interesting to remark that particles with masses and widths very similar to those 

above for the a and the /o(980) were predicted |n| as part of a multiquark qqqq nonet on the 

basis of the MIT bag model. Hence, even though they do not give rise to formally leading 

7T7T amplitudes in the — scheme, the picture has a good deal of plausibility from a polology 

1 

point of view. It is not hard to imagine that some — — subleading effects might be important 
at low energies where the QCD coupling constant is strongest. 

Other than requiring \R®\ < — we have not attempted to fit the puzzling experimen- 



tal results in the the /o(980) region. Recent interesting discussions are given in refs [18]. 
It appears that the opening of the KK channel plays an important role and furthermore, 
additional resonances may be needed. In this paper we have restricted attention to the 
7T7T channel (although the effective Lagrangian was written down for the case of three light 
quarks). Clearly, it would be interesting to study the /o(980) region in the future, according 
to the present scheme. 
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4.1 Imaginary part and Phase Shift 

Finally, let us discuss the imaginary piece I®. In the leading — limit the imaginary part 
vanishes away from the singularities at the poles, whereas Rq has support all over. This 

suggests that we determine an approximation to I® from the — leading Rq using dispersion 

■N't c 

theory, rather than getting it directly from the tree amplitude with the regularization of the 

form ( [2.21 ). The latter procedure picks up pion loop contribution to the p — propagator, for 

example, but misses very important direct pion loop contributions. A dispersion approach 

will include both. In the low energy region, we can proceed more simply by just using the 

unitarity formula ( |2.5| ) directly. Up until the KK threshold it seems to be reasonable to 

approximate the elasticity function t)q(s) by unity [18]. Strictly speaking tjq(s) may depart 

from unity at the 47r° threshold^] of 540 MeV. In Fig. 12 we show I® obtained from ( |2.5D 

on the assumption t]q = 1 for several values of G' . Both signs in front of the square root 

are displayed. Of course, the correct curve should start from zero at threshold (— sign in 

front of the square root). Continuity of I® (s) would at first appear to suggest that we follow 

along the lower curve. In order to go continuously to the upper curve it is necessary that the 

argument of the square root vanish at some value of \fs. With the approximation T]q = 1, 

1 

this vanishing occurs if \Rq\ is exactly — . In Fig. 12, the discontinuity in the \/s = 540 MeV 

2 1 
region is extremely sensitive to tiny departures of \Rq\ from — . However, both experiment 

d.5 1 

and the expectation that 1 > [] suggest that beyond yj~s 540 MeV we should actually 

ay s 

go to the upper curve (+ sign in front of the square root). This can be accomplished without 

6 It is amusing to note that each of the low energy resonances, i.e. <r(530) and /o(980), are located just 
below threshold; for the a it is the 4n threshold while for the /o(980) it is the KK threshold. 

It/ 

7 In potential theory Wigner j2o) has shown that — l — > — — where a is the approximate interaction radius 

(X-y S fj 

and (3 is the pion velocity in the center of mass. Strictly speaking, for a < 1.7 fm the lower curve is also 
allowed. 
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violating continuity by assuming that t/q is not precisely one. For the curves shown all 
that is required is a decrease in t]q of not more than 0.04. Alternatively, we could choose 
parameters so that R® reaches 0.5 precisely. Then the fit at higher energies is slightly worse. 
The corresponding three curves for the phase shifts are shown in Fig. 13. The discontinuity 
should be smoothed over in accordance with our discussion above. The agreement with 
experiment is quite reasonable up to about 860 MeV. We did not go beyond this point for 
the purpose of obtaining the phase shift because we are neglecting the KK channel which 
becomes relevant in the computation of the imaginary part. 

5 Summary and discussion 

In the leading large N c approximation to QCD, irn scattering corresponds to the sum of an 
infinite number of tree diagrams which can be of the contact type or can involve resonance 
exchange. This can only be a practically useful approximation if it is possible to retain just 
a reasonably small number of terms. The most natural way to do so is, of course, to consider 
contact terms with as few derivatives as possible and exchange terms with resonances having 
masses less than the extent of the energy region we wish to describe. In this paper we have 
carried out an initial exploration of this program in a step by step way. The first step is to 
include only the well known chiral contact term which reasonably describes the scattering 
lengths. However this amplitude badly violates partial wave unitarity bounds (seen most 
readily in the I = L = channel, see Fig. 1) at energies beyond 500 MeV. We observed 
that the introduction of the p meson dramatically improved the situation, delaying drastic 
violation of the unitarity bound till around 2 GeV (see Fig. 2). We noted that this effect 
could be nicely understood as the result of an extra contact term which must be present 
when the p is introduced in a chirally invariant manner. Furthermore, this feature holds in 
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the strict large iV c limit, i.e., without including the phenomenological regularization (2.21). 
The observed cancellation encouraged us to investigate the possibility of a more general local 
cancellation, due to inclusion of all (large iV c leading) resonances in the energy range of 
interest and, perhaps, higher derivative contact terms. The program is sketched at the end 
of section 2. 

Taking the large N c approach as well as the standard qq spectrum literally, we argued 
that the next group of resonances whose exchange contributes to the leading amplitude should 
comprise the /o(1300), the /2(1270) and the ^(1450). We observed (section 3) that there was a 
tendency for these to cancel among themselves; for example the crossed-channel exchanges of 
the /?(1450) tended to cancel against those of the /2(1270). In our analysis, the complications 
due to enforcing chiral symmetry and using the full spin 2 propagator were taken into account. 
However, the cancellation with both the ir + p and next group was not sufficient to satisfy the 
unitarity bound \Rq\ < — in the energy range till 1.3 GeV. An allowed leading N c way out - 
by adding four derivative contact terms - was thus investigated. This enabled us to restore 
unitarity till about 1.0 GeV (see Fig. 7). The drawback was that Rq(s) dropped off rather 
precipitously afterwards, which would make a local cancellation scheme very unstable. 

As a more physically motivated alternative we investigated, in section 4, the possibility 
of including scalar resonances having masses less than 1 GeV. These are presumably not of 
the simple qq type and hence their exchange should be of sub-leading order in the large N c 



limit. An interesting interpretation gives these particles a qqqq quark structure [19]. Then a 
somewhat narrow state like the /o(980) is expected together with a very low mass and very 
broad state like the old a meson. (Both should belong to a 3— flavor nonet). It was found 
that the /o(980) particle did not help much in restoring unitarity. In the experimentally 
puzzling region close to 980 MeV it is, however, expected to play an extremely important 
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role. On the othe hand, the further introduction of the other scalar, which we denoted as 
the <r(530), treated with a phenomenological regularization parameter (see ( |4.1| )) enabled 
us to satisfy the unitarity bound all the way up to 1.3 GeV (see Fig. 10). Thus, if low 
energy scalars are included, the proposed local cancellation may be a viable possibility. The 
imaginary part of the partial wave amplitude, Iq(s) was also computed from the unitarity 
relation ( |2.5[) and found to lead to a phase shift <5§ in reasonable agreement with experiment 
until about 860 MeV. Beyond this point, the effect of the opening of the KK channel must 
be specifically included. 

There are many directions for further work. 

i. The most straightforward is the investigation of different channels. For example, con- 
sidering 7T7T — * KK and KK — ► KK should enable us to study the interesting KK 
threshold region in a more detailed way. Looking at channels which don't communicate 
with 7T7T would enable one to focus on particular resonance exchanges. 

ii. The greatly increasing density of levels as one goes up in energy clearly indicates that 
there is a limit to how far one can go with the kind of microscopic approach presented 
here. It is expected that at energies not too much higher than the 1.3 GeV region 
this analysis should merge with some kind of string- like picture |§. In that region 
the question of the validity of the — expansion and a possible local cancellation can 
presumably be approached in a more analytical manner and interesting models can be 
studied |21]. Here we have tried to follow a phenomenologically oriented path, assuming 
only chiral dynamics in addition to the — framework. 



Hi. One can also imagine a kind of Wilsonian effective action [22] with which the present 



approach can be further discussed. This should allow the systematic calculation of 
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loops but would be extremely complicated in practice. 
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Appendix A 

Here we list the kinematic conventions for tttt scattering. The invariant amplitude for 7Tj + 
TTj — > vrfc + 7Ti is decomposed as: 

SijS k iA(s, t, u) + S ik 5jiA(t, s, u) + SuS jk A(u, t, s), (A.l) 

where s,t and u are the usual Mandelstam variables obeying s + t + u = 4m^.. Physical values 
lie in the region s>4m^, t<0,u<0. (Note that the phase of ( |A.1| ) corresponds to simply 
taking the matrix element of the Lagrangian density of a four point contact interaction). 
Projecting out amplitudes of definite isospin yields: 

T°(s,t,u) = 3A(s,t,u) + A(t,s,u) + A(u,t,s), 

T l (s,t,u) = A{t,s,u)-A(u,t,s), 

T 2 (s,t,u) = A{t,s,u)+A(u,t,s). (A.2) 
In the center of mass frame: 



2 , 2\ 
-I* ), 



s = A(p w + m„ 
t = -2pl(l-cos6), 

u = -2j4(1 + cos9), (A.3) 

where p n is the spatial momentum and 9 is the scattering angle. We then define the partial 
wave isospin amplitudes according to the following formula: 



If (a) = (l"4^) j^dcoaOP l {cose)T I (s,t,u). (A.4) 
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the Roy curves for R®. 
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4 The solid line is the ir + p contribution to the 1 = 2, L = real part. The dot-dashed 
line is the it + p contribution to the / = 0, L = 2 real part. 
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6 Sum of all contributions in Fig. 5. 
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Solid line: vr + p(770) + / (980) + / 2 (1275) + /o(1300) + p(1450). Dashed line: without 
/o(980). 
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Fig. 11 The low energy structure for current- algebra, solid-line; ir + p, dashed-line; everything, 
dotted-line. 
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Fig. 12 Estimated imaginary part Iq. Upper curves + sign in front of the square root in fl2.5|), 
lower curves, - sign. Solid line: G' = 380 MeV. Dashed line: G' = 440 MeV. Dotted 
line: G' = 470 MeV. 
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Fig. 13 Estimated phase shift <5§. Solid line: O = 380 MeV . Dashed line: G' = 440 MeV. 
Dotted line: G' = 470 MeV. 
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